arXiv: 1505.07678v2 [nucl-th] 3 Sep 2015 


Microscopic evaluation of the hypernuclear chart with A-hyperons 

E. Khan,' J. Margueron,^ F. Gulminelli,^ and Ad. R. Raduta^ 

^Institut de Physique Nucleaire, Universite Paris-Sud, IN2P3-CNRS, F-91406 Orsay cedex, France 
^Institut de Physique NucUaire de Lyon, University Claude Bernard Lyon 1, 

IN2P3-CNRS, F-69622 Villeurbanne cedex, France 
^ CNRS/ENSICAEN/LPC/Universite de Caen Basse Normandy, UMR6534, F-14050 Caen c^dex, France 
^IFIN-HH, Bucharest-Magurele, POB-MG6, Romania 
(Dated: September 4, 2015) 

Background: A large number of hypemuclei, where a considerable fraction of nucleons is replaced by strange baryons, and 
even pure hyperonic species are expected to be bound. Though, the hypernuclear landscape remains largely unknown because 
of scarce constraints on the NY and YY interactions. 

Purpose: We want to estimate the number of potentially bound hypemuclei. In order to evaluate realistic error bars within the 
theoretical uncertainties associated to the spherical mean-field approach, and the present information from already synthetized 
hypemuclei on the N — Y and Y — Y channels, we limit ourselves to purely A-hypemuclei, to magic numbers of A’s (for Z < 
120 and A <70), and to even-even-even systems. 

Method: We consider a density functional approach adjusted to microscopic Bmckner-Hartree-Fock calculations, where the 
AA term is corrected in a phenomenological way, to reproduce present experimental constraints. Different models which 
strongly deviate at large densities, but giving the same bond energy, are generated in order to take into account the uncertainties 
related to the high density equation of state. 

Results: The number of bound even-even-even A-hypemuclei is estimated to 491680 ± 34400. This relatively low uncertainty 
is due to the fact that the well constrained low density and highly unconstrained high density behavior of the energy functional 
turn out to be largely decoupled. Results in A-hypemuclei appear to be almost independent of the choice for the high-density 
part of the AA interaction. The location of the A-hyperdriplines is also evaluated. Significant deviations from Iron-Nickel 
elements can be found for A-hypemuclei with the largest binding energy per baryon. Proton, neutron and A-hyperon magicity 
evolution and triple magic A-hypemuclei are studied. Possible bubbles and haloes effect in A-hypemuclei are also discussed. 

Conclusions: The present results provide a first microscopic evaluation of the A-hypemuclear landscape. They shall benefit 
from the more and more accurate design of the A-based functionals. The measurements of A and multi-A hypemuclei together 
with additional constraints of the SA and SS interaction is mandatory to improve such critical information. 

PACS numbers: 21.80.+a, 21.60.Jz, 21.10.Dr 


I. INTRODUCTION 

The study of hypernuclei benefited from a great scientific 
interest since the ’60s Gna. Various hyperon-nucleus interac¬ 
tion El, Skyrme Hartree-Fock iHni, relativistic mean-field 
lfT2UT9l . generalized liquid drop ll20l - E4ll and, G-matrix ll25l - 
Eil models have been considered in order to underpin the 
hyperon-nucleon interaction from hypemuclei, to test the ex¬ 
istence of bound hypemuclei, the stability of nucleonic cores 
against hyperon addition, or the occurrence of exotic pure 
strange systems or halo structures. With the advent of new 
dedicated experimental programs such as J-Parc in Japan or 
PANDA at FAIR, the study of the hypemuclei structure en¬ 
joys a revived interest Il28l - l30l . 

Present facilities can only produce single and double A 
hypemuclei in a limited domain of mass. However, a gen¬ 
eral understanding of the specificity of hypernuclear structure 
with respect to nuclear structure requires the evaluation of the 
global hypernuclear chart, with strangeness as the third di¬ 
mension cma]. Some recent works 1291 address this 
problem within phenomenological mean field models. How¬ 
ever, the uncertainties on the hypernuclear chart associated to 
the choice of the functional in the strangeness sector are dif¬ 
ficult to evaluate, since present experimental data on hyper¬ 
nuclei are scarce, and therefore a large arbitrariness is associ¬ 


ated to the modeling of the hyperonic energy functional. To 
minimize such uncertainties, it is important to use as much 
as possible microscopically founded energy functionals from 
Brueckner or Dirac-Braeckner calculations as well as the few 
available experimental data. The present work is an attempt 
towards that direction. 

The regular nuclear landscape, including an estimation of 
the uncertainties on its limits defined by the drip-lines, has 
been only recently microscopically studied E21 j owing to the 
growth of calculation capacities. This is also due to the re¬ 
cent use of evaluation methods mM for the uncertainties 
generated by the only input of such microscopic calculations: 
the nucleon-nucleon (NN) energy functional. The resulting 
uncertainty on the number of bound nuclei is typically of the 
order of 7% l32l . 

The present work aims to generalize such a study by de¬ 
signing the limits of the A-hypernuclear landscape, and eval¬ 
uating the uncertainties associated with the limited empirical 
information on the hyperonic functionals. For this purpose 
microscopic calculations of nuclei and A-hypernuclei are per¬ 
formed using the energy density functional approach. We use 
a Skyrme functional optimized on nuclear physics data for the 
NN channel, while AN and AA functionals fitted on micro¬ 
scopic Brueckner-Hartree-Fock (BHF) calculations are used, 
with modifications in order to take into account A-hypernuclei 
data constraints. 
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When dealing with multi-strange system, a very important 
issue is the possible presence of hyperons other than A. In par¬ 
ticular, it was recognized since the early nineties |[I61IIS|351 
that the A + A^Z+N decaying channel could play a major 
role in multi-strange systems, because of the attractive charac¬ 
ter of the E potential in nuclear matter, which could lead to the 
appearance of E hyperons already for a strangeness number 
|5| > 8. The present most accepted value for the E potential 
in symmetric nuclear matter Vs = —14 MeV ll^ is much less 
attractive than the value proposed in those earlier works m, 
which pushes the A number threshold for E contribution far¬ 
ther in strangeness. Moreover, it is clear that the contribution 
of E will crucially depend on the E-E and E-A interactions, 
which are presently completely unconstrained. The inclusion 
or removal of the Fock term in relativistic approaches seams 
to have, as well, an effect on the appearance the various kinds 
of hyperons in uniform matter BTI . Because of these uncer¬ 
tainties, we have chosen to limit ourselves to A- hypernuclei 
in this work, since A’s are the only hyperons to be (relatively) 
constrained by experimental data. 

Section II details the energy density functional used, fo¬ 
cusing on the way to design the AA component. Section III 
is devoted to the determination of the corresponding parame¬ 
ters of the AA functional. The Hartree-Fock calculations for 
A-hypernuclei are performed in Section IV. Since in the hy- 
peronic channel the spin-orbit interaction is expected to be 
very weak 1^ . the A-hypernuclear charts, location of the 
A-hyperdriplines and the estimation of the number of bound 
even-even-even A-hypernuclei are evaluated for hypernuclear 
number A=0, 2, 8, 20, 40 and 70. In Section V, the gross 
properties of A-hypernuclear structure are analyzed, namely 
the evolution of the energy per baryon as a function of the 
AN and AA functionals. A study of magicity evolution in 
A-hypernuclei, as well as of possible bubbles and haloes ef¬ 
fects is also undertaken. Finally, Appendix A details the link 
between the bond energy and the AA functional whereas Ap¬ 
pendix B provides an update of the strangeness analog reso¬ 
nances in multihyperons A-hypernuclei. 


II. DENSITY FUNCTIONAL THEORY FOR 
A-HYPER-NUCLEAR MATTER AND A-HYPER-NUCLEI 

We consider a non-relativistic system composed of interact¬ 
ing nucleons N and lambdas A. The total Hamiltonian reads, 

H — tm+ t/^ ++ H/^ + (1) 

where Ta is the kinetic energy operator and Hab the interac¬ 
tion operator acting between A and B {=N and A). We work 
in the mean-field approximation where the ground state of the 
system is given by the tensor product, |<I> 7 v) ® I'I’a); where 
|<I>a) = n,a^|—) is a Slater determinant of nucleon states and 
IOa) = Ilxa^\—) is a Slater determinant of lambda states. 
The total Hamiltonian Q can be expressed in terms of the 
nucleons (/) and lambda (A) states as, 

«=Et+Efi+5E'=!;''+ { E'Ta+E'l''- m 

i A ^ ij X.fi X,i 


In the following, we will consider the density functional 
theory which allows relating in a direct way the BHF predic¬ 
tions for uniform matter to the properties of hyper-nuclei. 

A. Energy-density functional deduced from BHF 

In the present study of A-hypernuclei and nuclear matter we 
use a density functional which has been determined directly 
from BHF theory including nucleons and A-hyperons ll2^IZ7l . 
The total energy density £(pa,Pa) is related to the energy 
per particle calculated within the BHF framework, cbhf, as 
e(PA',PA) = (Pn + PA)eBHF{PN,PA) and is decomposed in 
different terms, 

. , tP- fP 

£\Pn,Pa) = ^ —Tat-I--—' fA + ewlPiv) 

Zlfl\ 

+£na{Pn,Pa)+£aa{Pa), (3) 

where, in infinite nuclear matter, the kinetic energy densities 
Tat and Ta are simple functions of the matter density: T,- = 
p^i^^, with gi = 4(2) for i = N{A). 

In the nucleon sector, we use the SLy5 parametrization 
of the phenomenological Skyrme functional including non¬ 
local and spin-orbit terms, since it can correctly reproduce 
the properties of stable and exotic nuclei lf39l . In the case 
of the strangeness sector, the spin-orbit interaction is known 
to be small lf38l . and is therefore neglected. The local den¬ 
sity dependence of the NA component of the energy den¬ 
sity, £na{PNjPa), is solely adjusted to the BHF predictions. 
To pin down the NA coupling, the following energy den¬ 
sity is defined: (pa- + Pa)^^//^(Pa,Pa) — PNeBHFiPN,^) — 
PAeBHF{0,PA), where eBHF{PN,PA) is the BHF energy per 
baryon in an infinite hyper-nuclear matter calculation. It is 
parameterized in terms of the nucleon and hyperon densities 
as 12611221, 

£na{Pn,Pa) = -f\ (Pa')Pa'Pa +/2(Pw)PaPa^^) (4) 

where the first term physically corresponds to the attractive 
NA interaction, corrected by the presence of the medium 
given by the function fi, and the second term is induced by the 
repulsive momentum dependent term of the A potential (con¬ 
sidering the low-momentum quadratic approximation), also 
corrected by the medium through the function fx- In the pres¬ 
ence of the attractive AA interaction, the term Caa is solely 


determined by the hyperon density as EH, 

£AA{pA) = -h{pA)pi- (5) 

The functions /,■ are given by the polynomial forms, 

/i (pa) = «! - ajpN + a-ipl, (6) 

fiiPN) = a4 - asPN + OePA- (7) 

/3 (Pa) = a? - 0!8 Pa + a9PA- (8) 


In principle, the functions f\ and fi depend on the densities 
Pa and Pa associated to conserved charges in the medium. 
However, since nucleons are the dominant species, even in 
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TABLE I. Parameters of the f, functions, see Eqs. for the functionals DE-NSC89, DE-NSC97a, DF-NSC97f. 


Force 

«! 

«2 

as 


as 

CC6 

a^ 

ag 

09 

DF-NSC89 12611271 

327 

1159 

1163 

335 

1102 

1660 

0 

0 

0 

DF-NSC97a j27l 

423 

1899 

3795 

577 

4017 

11061 

38 

186 

22 

DF-NSC97f (27) 

384 

1473 

1933 

635 

1829 

4100 

50 

545 

981 


the presence of lambdas, the dependence on p\ is neglected 
in these functions. In the case of it trivially depends on Pa 
only, because it impacts the part of the functional referring to 
pure lambda matter. 

Different NA potentials have been parameterized which 
equally well fit the scarce NA phase shifts, see for in¬ 
stance discussion and references in Ref. ioi. The present 
study is based on three of them for which a density func¬ 
tional has been derived, based on BHF predictions ll^ IZTl . 
namely DF-NSC89, DF-NSC97a and DF-NSC97f. The func¬ 
tional DF-NSC89 is based on the Nijmegen soft-core poten¬ 
tial NSC89 mi, while the functionals DF-NSC97a and DF- 
NSC97f are based on two of a series of six different hyperon- 
nucleon potentials which equally well reproduce the measured 
scattering lengths in the AN and T.N channels, and correspond 
to different values of the magnetic vector meson ratio a which 
cannot be constrained by the phase shift information Il42ll43l . 

Although more recent functionals have been parameter¬ 
ized M, by adding for instance isospin degree of free¬ 
dom, the functionals considered in this study already repre¬ 
sent a good sample of the uncertainty generated by the lack 
of empirical information in the strangeness sector. This af¬ 
fects microscopic approaches like BFIF |[8l|26][27l, although 
to a lesser extent than fully phenomenological mean field 
models El III HHIH HMH125] ElHlol. Indeed the models 
DF-NSC97a and DF-NSC97f correspond to the two extreme 
choices for the unconstrained a parameter Il42ll4^ . leading to 
the softest and stiffest equation of state, respectively. It should 
be noted that no experimental information is available on A-A 
scattering, meaning that these phenomenological bare inter¬ 
actions are completely unconstrained in the A-A channel. For 
this reason, NSC89 does not contain any A-A interaction. The 
NSC97 models assume for this channel a simple SU(3) ex¬ 
tension of the original Nijmegen potential models to multiple 
strangeness S = —2. The values for the parameters ai — Og of 
the functions fi — fj, are given in Table |I] for the functionals 
DF-NSC89, DF-NSC976a and DF-NSC97f 


The single particle energies in uniform matter are deduced 


from the energy functional (31, as = 


unif . 


unif 






, where 


y7^iPN,PA) = v7 + 


sky 


deNA 
dpN ’ 


(9) 


being deduced from the Skyrme functional Il44ll45l . and 


vT^iPN.PA) 


^Ena ^£aa 
dpA dpA 


( 10 ) 


In the following, the hyperon potential is decomposed 



FIG. 1. (Color online) Potential {pm,Pa) a function of 

the nucleon density p^ (in units of the saturation density po) for the 
functionals DF-NSC89, DF-NSC97a and DF-NSC97f without A (a) 
and with Pa = 0.03 fm^^ (b). 


into two terms, using Eqs. and 

0, 


^{N),umf ^ ^ {Pn)Pn + :^f2{PN)PNpf^ ( 11 ) 

^{A),umf ^ ^ -laiPA + 3aspl - Tagpl- (12) 

a Pa 


The term stands for the contribution of the nucle¬ 
ons to the hyperon potential, while the term repre¬ 

sents the direct contribution of the hyperons to their own po¬ 
tential. The properties of the potential y^l’""'^ are analyzed 
in Fig. [^imposing p„ = Pp. On the left panel, the poten¬ 
tial v' 


(Npunif . 


is displayed without A particles, while on the 
right panel a small amount of A is considered, correspond¬ 
ing to a representative average A-density in single A hypernu¬ 
clei (see Sec. IIIAi. It is expected, from experimental single 
A-hypernuclei data, that the potential v7 is about -30 MeV 


at saturation density 112611 . This empirical condition is sat¬ 
isfied for the three functionals DF-NSC89, DF-NSC97a and 
DF-NSC97f at saturation density, as shown on the left panel 
of Fig. [T] For a fixed and small amount of A, the potential 
v^)’“"‘Ms attractive at large densities for all the functionals. 


















4 


due to the term in the function f\. Already in the late 
eighties it was clear that some repulsion is necessary at high 
densities to explain hypernuclear data 0. For a finite amount 
of A, the as term in function /2 gives a repulsive contribution, 


and can compensate the attractive a$ term if Pa > ^ 


5 CCS J 


For the functionals DF-NSC89, DF-NSC97a and DF-NSC97f, 
this occurs for Pa > 0.27, 0.09 and 0.15 fm^^, respectively. 
These substantially different numbers reveal the large uncer¬ 
tainties at high density in the predictions of the considered 
functionals. These three functionals have been used before 
in Skyrme-Hartree-Fock calculations of single A-hypernuclei 
from C to Pb, and the resulting A-single particle levels have 
been confronted in detail to experimental data IIZTll . The sin¬ 
gle particle levels of the heaviest A-hypernuclei, which should 
in principle be reasonably well described in a mean-field cal¬ 
culation, were seen to be systematically over-bound by DF- 
NSC97a, and slightly under-bound by DF-NSC89, while DF- 
NSC97f nicely reproduces the data. In this sense, we can con¬ 


sider that the A potentials in uniform matter presented 

in Fig.l span the range compatible with A-hypernuclear data, 
even if the functionals appear to be probably too attractive at 
high density lIZTl . 


These functionals have also been extended to other hyper- 
onic channels, namely T, and E, assuming SU{3) symmetry, 
and exploited in BFIF matter calculations ll46ll . However, 
these extrapolations of the N — A channel to the rest of the 
octet, have been severely questioned in recent years due to 
the difficulty of standard functionals to predict the existence 
of very massive neutron stars MM- For this reason we 
will not consider these extensions here. It is however in¬ 
teresting to note that all functionals predict that, for a suffi¬ 
ciently high number of A, the A +A —> E -|-A decay channel 
should open and E should contribute to the multistrange nu¬ 
clei. The threshold A fraction for E appearance is estimated to 
be Pa/Pn = 0.17 for DF-NSC89 the only functional that pre¬ 
dicts a E potential in symmetric matter Ve = —15 MeV lIZTl . 
compatible with experimental data 1^ . It should be noted 
that this criterion is largely influenced by the value of the E 
potential in symmetric matter, as expected, since for instance 
the threshold A fraction for E appearance is estimated to be 
p^/p^ = 0.15 and 0.08 fm^^ for the DF-NSC97a and DF- 
NSC97f functionals having more attractive E potentials M- 


This means that the results we will get for highly strange 
A-hypernuclei have to be taken with caution. In particu¬ 
lar, driplines associated to a lambda fraction Pa/Pn = 0.17 
have to be considered as a lower bound, because even more 
strangeness would be compatible with bound systems if E 
would be accounted for. 


The two terms and contributing to the po¬ 

tential (Eq. (10 1 ) are compared in Fig. as a function 
of the nucleonic density p^ and for different proportions of 
A. As expected, for a sufficiently large amount of A, the term 
^{N),unif repulsive at high nucleonic density. At low 

density the functionals DF-NSC97a and DF-NSC97f predict 


very close potentials but as the nucleonic density in¬ 

creases the predictions of these two functionals increasingly 


Pa=0.1 Ps Pa=0-2 Pn 




0 1 2 3 40 1 2 3 4 

Pk / Po Pn / Po 


FIG. 2. (Color online) Comparison of y^^’"”'-^ (top panels: a, b) 

and (bottom panels: c, d), as a function of the nucleon den¬ 

sity Pa (in units of the saturation density po) for the functionals DF- 
NSC89, DF-NSC97a and DF-NSC97f and a A fraction set to 10% 
(left panels: a, c) and 20% (right panels: b, d). 


deviate. Since the dominant term in Og is attractive, the po¬ 
tential yj^^’“"'^ will finally curve down at very high densities, 
but this occurs in physical situations which will be hardly met 
in nature. Indeed, for the functional DF-NSC97a, the poten¬ 
tial continuously increases up to p^ « lOpo and for a fraction 
of hyperons < 50%. For DF-NSC97f, the effect of Og is al¬ 
ready curving down the potential for p^ « 8po and 20% of 
A. In pure A matter, the potential becomes repulsive 

between po/2 and po. The functional DF-NSC97f predicts 
also a decreasing potential starting from 1.5po while this be¬ 
havior is pushed to much larger values (beyond 20po) with 
DF-NSC97a. 

These differences between the functionals clearly show that 
the high density properties of A-hyper-matter are largely un¬ 
constrained. From the phenomenological point of view, the 
existence of very massive neutron stars (about 2 M 0 mi) 
requires a strongly repulsive potential at high density, 
which corresponds to a low value of the a parameter, in terms 
of the elementary Nijmegen interactions. In principle, the 
very existence of such massive stars could be used as an ex¬ 
tra constraint for the functional, and was often considered in 
the relativistic mean-field literature (see for instance the recent 
work im and Refs, therein). However, due to the presence of 
different hyperon species in the core of massive neutron stars, 
with highly unknown interaction couplings, it is not straight¬ 
forward to convert this qualitative statement into a sharp con¬ 
straint on the potential. More specifically, it should be 
noted that no BHF calculation is presently able to reproduce 
the empirical observation of two solar mass neutrons stars. 
This contradiction between our knowledge of supra-nuclear 
matter and observation is usually called the hyperonization 
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puzzle ll50l . 

Finally, it can be noted from Fig. 2 that the term is 

much smaller than v^^’™‘^contributes to less than 

10% to the total potential In the functionals DF-NSC97a 

and DF-NSC97f, the contribution of the term in¬ 

duced by the AA interaction, to the properties of hyperonic 
matter is therefore expected to be rather weak. As discussed 
above, this result might not be entirely physical, since the A-A 
interaction in the NSC97 models is not fitted on experimental 
data, but only extrapolated from the N — A interaction. The 
AA term of the functional can thus be phenomenologically 
adjusted as explained in the next section. 


B. Phenomenological correction to the term e^A 


iii) Finally, we impose a relation between the AA interac¬ 
tion and the bond energy in ®He, obtained from the local 
density approximation, that we discuss hereafter ill. 

Condition i) imposes (% < 0. For convenience, we set 
Ctg = 0, giving minimal repulsion. Fig. shows that the con¬ 
tribution of the term appears at very high densities which 
are certainly never reached in neutron stars. 

Condition ii) and Eq. ([^give 

2a7 = 3a8PA = 3a8xpo, (15) 

meaning that the parameter x controls the high density behav¬ 
ior of A matter; the larger x, the softer the EoS with hyperons. 

Condition iii) is related to the bond energy ABaa which is 
defined as ||2^|Z7]| . 

ABaa(A) = -£(^-22)+2£(1-'z)-£(1aZ), (16) 


The term Caa defined by Eq. Q is generated by the AA 
interaction: only the functionals DE-NSC97a and DE-NSC97f 
for which the AA interaction has been included in the BHE 
calculation, have a non-zero Eaa 12611221 . 

However, large uncertainties remain on the AA interaction, 
as discussed above. The functionals introduced in sectioning 
based on an SU(3) extrapolation of the N — A bare interac¬ 
tion to the S = —2 channel iHTl l43l . do not lead to a satis¬ 
factory description of the binding energy of double-A hyper¬ 
nuclei na, which is the only empirical information that we 
have on A — A couplings. EBia. Moreover, Eig. shows 
that, though the DE-NSC97a and DE-NSC97f functionals in 
the A — A channel are very similar at low density, they drasti¬ 
cally differ above saturation. In view of the already mentioned 
hyperonization puzzle, it is therefore reasonable to consider 
that an even larger uncertainty in the high density behavior 
has to be associated to the term £aa- 

In the following, we therefore propose to modify the Caa 
term taking into account phenomenological arguments to deal 
with the uncertainty on the microscopic BHE results. In prac¬ 
tice, the values of the aj,as, OCg parameters of the original 
functional 0 are replaced by different values cty, ag, ag 
defining the new Eaa term, 

eAA = -(«7-«8PA + (%pl)pl, (13) 

and the new potential 


where is a nucleus with no hyperon, a is a single-A 
hypernucleus, and is a double-A hypernucleus. It should 
be noted that, experimentally, ®He and related hypernuclei 
provide the most accurate value of a bond energy with ABaa — 
1 MeV Il52ll5^. 

The relation between the bond energy and the functional is 
derived in Appendix [A] providing 


ABaa(A) « -2 


£aa(Pa(A)) 

Pa(A) 


(17) 


Injecting Eq. (13 1 into Eq. (17 1 , the following relation is de¬ 
duced. 


2pa(A) a7-asPA(A) =ABaa(A). 


(18) 


Using Eq. (15 1 and introducing xa(A) = Pa (A) / po, the aver¬ 
age A-density in double-A hypernucleus ^Z, we obtain from 
Eq. ([g. 


3x-2xa(A) XA(A)pod8 =ABaa(A). 


(19) 


We can thus express the new parameters dj and 018 as a 
function of x, xa(A) and ABaa(A). from Eqs. (15 1 and ([T^ 


as. 


a? = 


XABaa(A) 


3x —2 xa(A) 


xa(A)Po 


( 20 ) 


-(A),u.(f ^ _ 2 a^p^ + 2dgpl - 4dgpl (14) 


The determination of these parameters is based on the follow¬ 
ing phenomenological prescription: 


i) Since the AA interaction is expected to be repulsive at 
high density (to support for instance the observed 2 M 0 
neutron stars), the coefficient dg is taken as < 0 . 


ii) A parameter x is introduced, which represents the A 
density (in units of the saturation density po) where the 
A potential in hyperonic matter (I 4 i changes 

its sign and becomes repulsive. 


6:8 


ABaa(A) 


3x —2 xa(A) 


xa(A)p 


2 ' 

0 


( 21 ) 


The following expressions for the energy density and the 
potential are deduced from Eqs. (13 i and ([T4ll: 


eAA = -ypi(3xpo-2pA), 


A 


= -3dgpA{xpo- Pa)- 


( 22 ) 

(23) 


Since the parameters dj and dg are linearly related to the 
bond energy, the energy density Eaa also scales with the bond 
energy. The relation between cty and dg and the parameters x 
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and x\ is more complicated. Let us consider the limit where 
X 3> x\. In this case, the parameters a? and a% reduce to the 
simpler expressions; 


2xa(A)po 3xxa{A)p^ 


(24) 


It should be noted that is independent of x while is 
inversely related to x. Therefore the parameter ctg is mostly 
related to the prescription ii), modifying the high density part 
of the mean held, while a-j is almost independent of x. More¬ 
over, at the limit x ^ xa, the A-potential at the density Pa (A) 
is given by vj^^(pA = xa(A)po) « and the potential 

energy-density is given by eAA(PA = ■*a(A)Po) « e^A, where 

= -ABaa(A) (25) 

eAA = -0-5ABAA(A)xA(A)po. (26) 

The parameter is only related to the bond energy and is 
independent of x and xa(A), while the parameter e^A depends 
only of the bond energy and xa(A), independently of the value 
of X. At the limit where x ^ xa, the potential and the energy 
density at xa(A) depends only of the bond energy and xa(A), 
and are independent of x. More generally, the following rela¬ 
tions hold; 


vf\pA=XK{A)po) _ IXA 

" 3 X ^\\x) ) 

£aa(Pa =aa(A)po) ^ /f—VI 


(27) 

(28) 


where the correction goes like xa(A)/x. 

Eqs. (27 1 and ( [28] l show that the scaled AA potential 

and the scaled A energy Eaa/^aa solely depend on 
X and Xa, and are independent of the bond energy. A com¬ 
plete representation of the parameter space can therefore be 
obtained by representing these scaled quantities as a function 
of the A-density, and largely varying only the two parameters 
X and Xa. These normalized potentials (respectively energies) 
are represented in Fig. (respectively |^. 

According to condition ii), the mean field potential changes 
its sign for the density Pa = xpo. This is well observed in 
Fig. [3 Therefore, the effect of increasing the parameter x 
makes the A-mean field potential and energy density softer. 
In a similar but weaker way. Figs. andshow that the effect 
of the parameter xa is to soften the mean field potential and 
the energy density as xa increases. 

The potential and energy density behavior, displayed in 
Figs|^and[^ show that this choice of parameters spans a wide 
range of qualitative behaviors of the energy functional. 


C. Application of the energy-density functionals to 
A-hypernuclei 

For the sake of completeness, we recall here the results ob¬ 
tained in ll26l IZTll . also used in the present work, describing 



FIG. 3. (Color online) Normalized potential (Pa)/|V'^'^^*| for 
several choices of x and xa, see the caption. 



FIG. 4. (Color online) Normalized energy density eAA(pA)/|£AAl 
for the same parameters as in Fig.[^ 


the implementation of the Hartree-Fock approach in the A- 
hypernuclear case. Nucleons and lambdas in the medium ac¬ 
quire an effective mass which is generated by the momentum 
dependence of the interaction. The description of hnite nuclei 
requires to disentangle the momentum dependent part of the 
in-medium potential (the correction to the masses) from the 
momentum independent one, hereafter called the local part. 
The local energy density is defined by; 

tP' , ^ Pp ^ 

+ + (29) 
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TABLE II. Parameters of the A-effective mass given by Eq. j34| l for 
the functionals considered in this paper. 


Force 

Pi 

Pi 

P3 

Pa 

DF-NSC89 I261I27I 

1 

1.83 

5.33 

6.07 

DF-NSC97a W\ 

0.98 

1.72 

3.18 

0 

DF-NSC97f (27) 

0.93 

2.19 

3.89 

0 


which can be recasted as; 


^NN — ~ 


^loc _ ^ 

%A“ ^NA — 


3pNh^ f 1 

10 V 8N 

3pA/t^ /f 1 

10 1 gA 



(31) 

(32) 


where is derived from the Skyrme functional Il44ll45l and 
EjVA is given by Eq. 0- We have already seen that the nucleon 
contribution to the lambda potential, e^A, is much bigger than 
the A contribution, Caa- For this reason, the contribution of 
the AA interaction to the A-effective mass can be considered 
as a small corrections, and it has been neglected in Ref llll. 
We can therefore write: 


e^(PA) = eAA(pA). (33) 

In Eqs. and (j3^ the local part of the energy density re¬ 
quires the knowledge of the effective masses m'^ and m\. The 
nucleon effective mass is given from the Skyrme interac¬ 
tion El HU. The effective mass of the A particles is mainly 
generated by the momentum dependence of the NA interac¬ 
tion, and it can be deduced from the BHE calculations. The 
result is expressed as a polynomial in the nucleonic density 
Pa as 12^ . 


m\{pN) 

'«A 


= Pi - P2Pn + PsPa - P4Pa- 


(34) 


The values for the parameters pi _4 for the functional consid¬ 
ered here are given in Table |I^ and the density dependence 
of the effective mass is shown in Eig. Since the effective 
mass is only necessary in the description of finite nuclei, we 
have limited the densities to values around saturation densi¬ 
ties in Eig. 1^ We can again observe, on the left panel, that 
the three functionals show a qualitatively similar behavior for 
the A effective mass at low density, with the stiffest model DE- 
NSC97f showing as expected the strongest momentum depen¬ 
dence. 

Given the functional form of and m* for the nucleons 
and the lambdas, Eq. Q can be rewritten in a form which ex¬ 
plicitly disentangle the terms coming from the local operator 
p and the non-local operator f = Vp V as. 


e(PiV)PA5TA,TA) 


2m*j^{pN) 2ml{pN) 


+eMv(PA') + 4 a(Pa':Pa) + eAA(PA)- (35) 


Minimizing the total energy, defined from the density func¬ 
tional ®, and using the Skyrme model for the nucleonic 




FIG. 5. (Color online) Ratio of the effective and bare A mass 
™a(P^)/'”a ( 4 ) and 1 (b) as a function of the nucleon 

density Pa (in units of the saturation density po) for the functionals 
DF-NSC89, DF-NSC97a and DF-NSC97f. 


part, we obtain the usual Schrodinger equation (i = N,A), 


-V- 




2m* (r) 


V + Vi{r) - iWi{r){V x (j) (p/,a(r) 


= -ei,a(Pi,air), 

where the nucleon potential Va is defined as. 


(36) 


FA(r) = vr^(r) + 


tWA 


dpM V'^Kpa) 


Ta s/3 


\ 2mA 5 2mA 
The hyperon potential Va is given by. 


Pa 


(37) 


VAir) = v‘‘f- 


ma 


,{Pn) 


-1 


2/3 

2mA 


(38) 


As in the uniform case, the hyperon potential Va is decom¬ 
posed into two terms, Va = where. 


^A W=''a ■ 

,,(A)/ ^ (A).unif 
^A ('■)=^A ■ 


mA 


Upa-) 


0 


(40) 


The factor (mA/m^(pA) — 1) in Eq. (391 is displayed on the 
right panel of Eig. showing that it is quite large around sat¬ 
uration density, between 0.2 and 0.5. The modification of the 
lambda potential in nuclei with respect to the uniform poten¬ 
tial can be quite important as scales approximately with the 
kinetic energy density of the lambdas. 
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FIG. 6. (Color online) Density profiles in isotopes with the addi¬ 
tion of one or two A hyperons. Black (red) lines: nucleon (A) density. 
Full (dashed) lines: one (two) A is added to the ^He nucleus. Three 
different functionals, from left to right DF-NSC89, DF-NSC97a, and 
DF-NSC97f are considered. 


III. DETERMINATION OE THE AA EUNCTIONAL 


In the following we will use the previous formalism to per¬ 
form calculations in finite nuclei, and to determine the param¬ 
eters of the AA functional. In the case of single-A-hypernuclei 
there is no AA interaction, and the corresponding terms in 
the functional shall therefore not be considered. We therefore 
set Caa = 0 for single-A-hypernuclei, while for double and 
many-A-hypernuclei, we consider the original expression 0 
for Eaa- 


A. The average A-density in hypernuclei 


According to Eq. <[13’ the bond energy and the average 
density in the double-A hyper nucleus are closely related. 
The bond energy can be experimentally constrained, but no 
experimental determination of the average A-density exists 
yet. Fig. displays the nucleon and A densities in ^He and 
^He for the functionals DF-NSC89, DF-NSC97a and DF- 
NSC97f, using SFy5 for the nucleons. The average single and 
double-A-densities (Fig. show a moderate dependence on 
the model. 

A definition of the average A-density from the A-density 
profile is required in order to properly define the parameter 
.rA(A) . Fig. 1^ clearly shows that the A-density profile does 
not have a flat behavior at the center of a small hypernuclear 
system; as a consequence, the standard deviation from the av¬ 
erage will be quite large. Table [HI] displays the calculated av¬ 
erage double-A densities and standard deviations for the same 
functionals as in Fig. using three different ways to estimate 


TABLE III. Calculations of average A-densities in ^He (in fm ^) 
for the functionals DF-NSC89, DF-NSC97a and DF-NSC97f. 



DF-NSC89 

DF-NSC97a 

DF-NSC97f 

10^(Pa)a 

2.25 

3.69 

2.31 

102 (Pa)v 

2.64 

3.48 

2.65 

W^{Pa)t 

2.51 

3.55 

2.53 

IO^cta 

1.97 

2.89 

1.85 

lO^o-v 

1.85 

2.74 

1.76 

lO^CTj- 

1.90 

2.79 

1.80 


the average A-density: 


/„ \ Jd^rpAPA 

(41) 

\ Jd^rpAPN 

{Pa)N = r ,3 „ 

J d^rpA 

(42) 

^ Jd^rpApT 

\Pa)t — .,3 ^ 

j d^r Pt 

(43) 


where pr = Pn + Pa- The standard deviation is defined as 
O’; = ^{P/Oi ~ (Pa)?’ "'ith i = A, N or T. The comparison 
between the different ways to extract the average A-density 
in Table [nl| shows a 15% deviation. As anticipated, the stan¬ 
dard deviation is very large, almost of the order of the average 
value. It is therefore difficult to properly define an average 
A-density in ^He. 

A general expression for the average A-density in hyper¬ 
nuclei was given in Ref. Et], supposing that the A-density 
scales with the nucleon density. The following general ex¬ 
pression lIZTl was therefore proposed: 

pf”""'(A)«A^, (44) 


where A=1 for single-A hypernuclei and A=2 for double-A 
hypernuclei. Fig.j^shows the ratio of the microscopically cal¬ 
culated average A-density {Pa)t, over the general expression 
(44i. The general expression is rather well satisfied for large 
A, but a systematic deviation is observed for small A. The fit 
of the deviation is also shown in Fig.|^and we obtain. 


{Pa)i 




0.957-2.54A- 


(45) 


This shows that Eq. (44i of Ref. lIZTl is not valid for A < 
50. It should be noted that we have used the expression {Pa)t, 
see Eq. as a reference to estimate the average A-density. 
Similar results are obtained with {Pa)n for large nuclei, but 
the definition (Pa)a gives a systematic increase by a factor of 
about 2 compared to the two other definitions, in large nuclei. 


B. Empirical determination of the term Eaa 


The relation between the ABaa(A) parameters and xa(A) 
is given by Eq. (17 1 for a given A-hypernucleus A. However, 
while the bond energy can be determined from experimental 
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FIG. 7. (Color online) Ratios of the A density (eq.l|43^) over the ap¬ 
proximate relation (eq. j44| () as a function of baryon number for A = 1 
(red) and A = 2 (green). A fit is also shown. Three different func¬ 
tionals, from left to right DF-NSC89, DF-NSC97a, and DF-NSC97f 
are considered. 





FIG. 8. (Color online) Relation between the parameter XA(^He) 
and the bond energy set to 5 MeV for He extracted from the 
HF calculation with the functional DF-NSC89 with x = 1/2 and DF- 
NSC89 with x = 4. 


mass measurements, the average A-density have never been 
measured. We have seen in the previous section that the aver¬ 
age A-density depends on the functional and that the standard 
deviation in ^He, the nucleus for which the best measure¬ 
ment of the bond energy exists ll5^l5Tl . is almost comparable 
with the average A-density. Two different methods are used 
to set the value of xa: 

EmpA: The value of xa/aaHc) is fixed to an average value 
(=1/6) in ^He, independently of the functional. 

EmpB; The value of xa is optimized to obtain a fixed bond en¬ 
ergy («1 MeV, consistent with Refs. Il52l l5^ or « 5 
MeV, as suggested in Ref Ell) in He for each func¬ 
tional. 


In prescription EmpA, the relation between the bond energy 
and the functional is based on the local density approximation 
(see Eq. (17 i and App.[A]for more details). It is therefore in¬ 
teresting to calculate the bond energy which is obtained from 
the microscopic HE calculations, in order to estimate the ac¬ 
curacy of the local density approximation. Doing so for DE- 
NSC89, DE-NSC97a and DE-NSC97f, and varying the pa¬ 
rameter X from 1/2 to 4, the difference between the bound en¬ 
ergy set to determine the parameters and the one calculated 
from HE calculation in ^He is less than 20%. A larger dif¬ 
ference is found between the parameter xa =1/6 and the aver¬ 
age A-density in ^He. We conclude that the the prescription 
EmpA, based on the local density approximation, cannot pro¬ 
vide accurate parameters. This is related the analysis of the 
A-density profile in ^He which is not flat enough to allow 
for the local density approximation (see Tab.[nl|. 

In the present approach, given by the prescription EmpB, 
the value of xa is not fixed a priori but it is varied and cor¬ 
related with the bond energy determined from the HE calcu¬ 
lation. In such a way, the parameter xa is treated as a varia¬ 


tional parameter allowing to fit the bond energy in He. The 
correlation between the parameter xa and the bond energy 
AB/^{A = 6)^^ is shown in Eig. for the functional DE- 
NSC89 with X = 1 /2 and DE-NSC89 with x = 4. The results 
for the functionals DE-NSC89, DE-NSC97a and DE-NSC97f 
using this empirical prescription (called EmpB hereafter) are 
given in Tab. IV obtained from the adjustment to the bond 
energy. Since a large arbitrariness is associated to the AA 
functional, the simplest polynomial form that allows for the 
needed repulsion at high density, and corresponds to a bond 
energy of 1 MeV in ®He and neighboring hypernuclei is cho¬ 
sen. We will therefore use the functionals of Table IIV I for the 
calculations of this work. Allowing for this large variation of 
the X parameter, the largely unknown behavior at supersatura¬ 
tion is therefore decoupled from the behavior at the very low 
densities implied for the hypemuclei. 


Einally, Eig. 9 compares the potentials obtained 

from the functmnals DE-NSC89, DE-NSC97a and DE- 
NSC97f associated to the empirical prescriptions EmpBl and 
EmpB2, where the bond energy is fixed to 1 MeV. Around 
saturation density, the different models shown in Eig. [^pre¬ 
dict negative values ranging from -1 to -2 MeV. At densities 
higher than normal nuclear saturation density, the prescription 
EmpBl makes the potential y^^’""'^ much stiffer than the pre¬ 
scription EmpB2. It is also interesting to compare the empiri¬ 
cal models shown in Eig. |^with the initial ones displayed on 
Eig. 0 the empirical prescription EmpB 1 produces the stiffest 
potential while EmpB2 the softest. Our empirical prescrip¬ 
tion increases the exploratory domain of variation for the AA 
potential, as well as includes the initial potential. 
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TABLE IV. Prescription EmpB. We present the adjustment of the parameter Xf^{A) to the bond energy (pa 1 or 5 MeV in '*He), the values of 
the parameters a® and af, and the ratio of the A-density to the saturation density in He. 


Pot. AN 

DF-NSC89 

DF-NSC89 

DF-NSC97a 

DF-NSC97a 

DF-NSC97f 

DF-NSC97f 

Pot. AA 

EmpBl 

EmpB 2 

EmpBl 

EmpB 2 

EmpBl 

EmpB 2 

X 

1/2 

4 

1/2 

4 

1/2 

4 




A6aa(4 

= 6) = 5 MeV 



a? 

150 

80.82 

150 

70.64 

160.48 

100.33 

as 

1250 

84.19 

1250 

73.59 

1337 

104.51 

ASaa(6)"^ 

5.05 

4.86 

5.17 

5.03 

4.84 

5.00 

xa 

0.125 

0.2 

0.125 

0.23 

0.115 

0.16 




ABaa(4 

= 6) = 1 MeV 



a^ 

36.05 

25.53 

39.46 

22.85 

49.24 

35.25 

as 

300.48 

26.60 

328.81 

23.81 

410.32 

36.72 

ABaa(6)"^ 

0.97 

1.1 

0.95 

1.04 

0.98 

1.05 

Xa 

0.1 

0.125 

0.09 

0.14 

0.07 

0.09 




Pn / Po 


FIG. 9. (Color online) Comparison of the potential ’™''^(piv,PA) 
obtained with the prescription EmpBl and EmpB2 (see legend) as a 
function of the nucleon density Pf^ (in units of the saturation density 
Po) for the functionals DF-NSC89, DF-NSC97a and DF-NSC97f. 
We have fixed AS^a = 1 MeV and the fraction of A to be 20%. On 
the left panel, we have fixed x= 1/2, while on the right, x = 4. 


C. Calculations in A-hypernnclei 


In the following, we consider the three functionals DF- 
NSC89, DF-NSC97a and DF-NSC97f corrected with the em¬ 
pirical prescriptions EmpB 1 and EmpB2 associated to a bond 
energy of 1 MeV. We therefore consider 6 NA + AA func¬ 


tionals, corresponding to Table IV together with SLy5 
for the NN functional. Eor each NA functional (DE-NSC89, 
DF-NSC97a and DE-NSC97f) the first one (e.g. DE-NSC89- 
EmpBl) corresponds to the lowest hyperon density in nuclear 


matter (x=l/2), whereas the second one (e.g. DF-NSC89- 
EmpB2) corresponds to the largest hyperon density (x=4). 

It should be noted that the present work is focused on the 
uncertainty generated by the A related functionals on the num¬ 
ber of bound systems. The uncertainties generated by the NN 
functional itself have already been studied on the nuclear chart 
and is typically of 7% on the number of bound nuclei ll^ . 

In the present work, the A-hypernuclear charts are calcu¬ 
lated for even-even-even A-hypernuclei. Since spherical sym¬ 
metry is imposed, only magic lambda number A-hypernuclei 
are considered. It should be noted that a fully microscopic 
deformed approach such as Ref. ||8] is the most accurate one 
to predict driplines. However such a task is numerically very 
demanding, and it is out of the scope of the present paper. 

As mentioned above, the spin-orbit interaction is known to 
be weak in the A channel lIMI . a factor 100 lower than in the 
nucleonic sector na, according to experimental data. The A 
magic numbers are therefore expected to be close to the har¬ 
monic oscillator ones: 2, 8, 20, 40, and 70. Calculations are 
performed for A-hypernuclear charts corresponding to these 
specific numbers of A. 


IV. THE A-HYPERDRIPLINES 


The A-hyperdripline has been studied in Ref. Il26l with a 
similar HE approach, showing that the maximum number of 
bound A in an hypernucleus is about 1/3 of the number of 
nucleons. We aim here to provide a more general study of hy¬ 
perdriplines, namely also showing the evolution of the proton 
and neutron driplines with the number of hyperons. It is im¬ 
portant to stress that driplines associated to a lambda fraction 
Pa/Pn ~ 0.17 have to be considered as a lower bound, be¬ 
cause even more strangeness would be compatible with bound 
systems if E, which are neglected in the present approach, 
would be accounted for. 

Figure displays the microscopically calculated A- 
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FIG. 10. (Color online) A-Hypemuclear charts for magic A numbers 
calculated with the DF-NSC89+EmpB 1 functional. 


FIG. 11. (Color online) Minimal (upper line) and maximal (lower 
line) neutron number value for hyper nuclei with A=20 using the 6 
functionals described in Table llVI 


hypernuclear charts for A=0, 2, 8, 20, 40 and 70 using the DF- 
NSC89+EmpBl functional. Adding A’s to a nucleus increases 
the binding energy for A < 40, but conversely decreases it for 
A=40 and 70. This is due to the balance between the attractive 
AN and AA functionals and the progressive energy hlling of 
the A states in the mean A potential. Fig. [^also shows that 
with a large number of lambdas, the corresponding hypemu- 
clear chart is shifted towards larger N, Z values. This effect 
is mainly due to the AN functional and is related to the maxi¬ 
mum fraction of 1/3 of hyperons with respect to nucleons, as 
mentioned above. For instance the exotic '^'^Th core becomes 
bound in the presence of 70 A. A similar stabilization of ex¬ 
otic nuclear core has been predicted in IfTSl . 

To get a more accurate estimation of the dripline. Fig. o 
displays the dripline neutron numbers for each Z value in the 
case of A=20 and for the 6 (AAh-AN) functionals that we con¬ 
sider. The dripline is here defined when the chemical poten¬ 
tial becomes positive. The results are rather similar among 
the various functionals. This shows that despite the very large 
differences among these functionals both in the NA and AA 
channels, the error bars are not so large for the drip-line de¬ 
termination, contrarily to what could be expected. This may 
be due to the narrower density range corresponding to finite 
hyper-nuclei compared to hyper-nuclear matter, since the 6 
functionals provide similar potentials for subsaturation densi¬ 
ties. 

The number of bound even-even-even A-hypernuclei found 
for A < 70 and Z < 120 are given in Table [V] The disper¬ 
sion due to the uncertainty on the A-related functional is rather 
small; the average total number of even-even-even hypemu- 
clei for A=0, 2, 8, 20, 40, and 70 is 9770 ± 429. Interpolating 
the number of even-even-even hypernuclei on the whole hy¬ 
pernuclear chart, provides 59538 ± 4020 hypemuclei with a 
larger relative systematic uncertainties, which can be inferred 


from the interpolation procedure. The interpolation is linearly 
performed between the calculated magic A hypercharts. 

We stress that this interpolation is only done with the pur¬ 
pose of estimating of the total number of bound systems. In 
order to have a detailed description of the non-magic hyper¬ 
nuclear chart, taking into account deformations would be es¬ 
sential 111. 

The corresponding uncertainty is calculated from the dis¬ 
persion of the total number of even-even-even hypernuclei ob¬ 
tained with the six (AAh-AN) functionals. It should also be 
noted that magic nuclei are specihc nuclei in the hypernuclear 
channel and therefore the interpolation procedure is not opti¬ 
mal. Considering the uncertainty on the bond energy (between 
1 MeV and 5 MeV as explained in the previous section), the 
number of estimated even-even-even nuclei is 61460 ± 4300. 
The total number of hypernuclei, considering the odd ones is 
therefore 491680 ± 34400. If the uncertainty of the NN func¬ 
tional 1321 is included (here in a decorrelated way), the total 
number of hypernuclei with A < 70 and Z < 120 is 491680 ± 
59000. 

The relative uncertainty is therefore of about 4% on the 
magic A hypernuclear charts. It is determined using the mean 
value and the corresponding standard deviation considering 
the 6 (ANh-AA) functionals. It should be noted that the DF- 
NSC97a functional is more attractive than the other function¬ 
als, which is the main contribution to the uncertainties on the 
number of bound hypernuclei. If the bond energy requirement 
is changed from 1 MeV to 5 MeV, the variation among the 6 
newly derived functionals is about 5 %. All in all, it is safe 
to consider an upper limit of 7% uncertainty due to the AA 
and the AN functionals on the magic A hypernuclear charts. 
In the case of interpolated values, the uncertainty is of about 
7% considering the 6 (ANh-AA) functionals, plus 5 % from 
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TABLE V. Number of bound even-even-even A-hypemuclei for A < 70 and Z < 120 



DF-NSC89 

DF-NSC89 

DF-NSC97a 

DF-NSC97a 

DF-NSC97f 

DF-NSC97a 

Number of nuclei 

-l-EmpB 1 

-l-EmpB2 

-l-EmpB 1 

-l-EmpB2 

-l-EmpB 1 

-l-EmpB2 

A=0 

1578 

1578 

1573 

1573 

1578 

1578 

A=2 

1628 

1640 

1617 

1619 

1634 

1621 

A=8 

1647 

1644 

1692 

1680 

1750 

1749 

A=20 

1650 

1681 

1696 

1724 

1675 

1683 

A=40 

1713 

1736 

1961 

1972 

1722 

1716 

A=70 

1162 

1237 

1746 

1886 

1127 

1152 

Total 

9378 

9516 

10285 

10454 

9486 

9499 

Total interpolated 

56140 

57287 

64459 

65587 

56646 

56809 


changing the bond energy requirement from 1 MeV to 5 MeV. 
These values should be compared to the relative uncertainty 
of about 7% when the NN Skyrme functional is changed ll32l 
on the nuclear chart. It shows that the uncertainty from the A- 
related functionals is not significantly larger than the one from 
the NN functional. This is due to the focused range of densi¬ 
ties (i.e. below the saturation one) relevant for hypernuclei, as 
mentioned above. 

Coming back to the sensitivity of the hyperdriplines to the 
A + A —>■ N + Z decay channel, we mention that within the 
generalized liquid-drop model proposed by Samanta 12^ and 
straightforward to use it comes out that the n-rich frontier is 
almost unaffected by considering Es in addition to As while 
the n-poor frontier is shifted to lower N-values. Quite re¬ 
markable, according to Ref. Il24l along the n-poor dripline 
strangeness consists almost entirely out of while along the 
n-rich dripline strangeness is to a large extent made of As. 
The broadening of the hypernuclear chart is proportional with 
the strangeness fraction. The same qualitative behavior is ex¬ 
pected to manifest also in the case of the present model. 


V. A-HYPERNUCLEAR STRUCTURE 
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FIG. 12. A-Hypemuclei with maximum binding energy per 
baryon as a function of A for the DF-NSC89-l-EmpB 1 (a), DF- 
NSC89-l-EmpB2 (b), DF-NSC97a-l-EmpBl (c), DF-NSC97f+EmpB 1 
(d) functionals 


A. Binding energy: fusion and fission of A-hypernnclei 

It may be relevant to study how the presence of hyper- 
ons impacts the most bound hypernuclei per baryon, which 
is known to be in the Fe-Ni region in the case of nuclei. 

Fig. [redisplays the evolution of the largest B/A value as a 
function of the number of hyperons in hypernuclei calculated 
with 4 (ANh-AA) functionals. A remarkable agreement be¬ 
tween the four functionals considered is observed, which may 
suggests that the theoretical uncertainty on the problem of hy¬ 
pernuclear binding is relatively under control, at least in the 
low density region corresponding to nuclei. It would be very 
interesting to know if this agreement is kept using relativistic 
functionals as in Ref. Il29ll . In the case of DF-NSC89H-EmpBl 
and for A=2 or 8, the optimal Z value is not much changed 
compared to the Nickel-Iron area. For a number of hyperons 
from 20 and larger the most bound hypemucleus is obtained 


for larger Z, typically for Zr, Ce and Pb. The results are qual¬ 
itatively similar with the other functionals, except in the case 
of DF-NSC97aH-EmpB 1 and DE-NSC97f+EmpB 1, for 2A hy- 
pemuclei: the addition of only 2A makes Si and, respectively, 
Ti the most bound nucleus in terms of binding energy per 
baryon. This result shows the non-negligible impact of hy¬ 
perons on the binding energy of the system. 


B. Magicity 

Several signals of the evolution of the magic gaps along the 
nuclear chart have been obtained these last decades EH. It 
may therefore be relevant to extend prediction of magicity to 
the hypernuclear chart. The two protons or two neutrons gap 
are known to be a relevant quantity in order to provide a first 
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FIG. 13. The 2 protons gap in the bound Nickel isotopes for A=0 
and A=40, defined as S 2 p = S 2 p{A,Z) — S 2 p{A + 2,Z + 2), calculated 
with the DF-NSC89+EmpBl functionnal 



FIG. 14. (Color online) Baryonic densities for the triply magic 
208+Apb hypemucleus with A=2,8,20,40,70 


insight on magicity features in nuclei. 


It is known that shell effects can vanish in the case of very 
neutron rich nuclei Il56ll57l because of the smoothness of the 
neutron skin triggering a weakening of the spin-orbit effect, 
and therefore its corresponding magic numbers. This is the 
case for hypernuclei with a low number of A. But in the case 
of A=40, Figure 13 shows a restoration of the Z=28 magic¬ 
ity for the very neutron-rich hypernickels with 70 < N < 80: 
the spin-orbit weakening of the Z=28 magic number for very 
neutron-rich hypernuclei is reduced by the presence of hyper- 
ons. A similar effect is observed for A=70 and the N=186 
magic number: it is restored for 82 < Z <98. The present re¬ 
sults were obtained with DF-NSC89-i-EmpBl and similar re¬ 
sults are obtained with the other functionals. A weakening of 
the N=28 shell closure is found for A < 8 in agreement with 
the previous results. Also a Z=40 weakening is observed for 
these hypemuclei. 


Fig. [^displays the neutron, proton and hyperon densi¬ 
ties for the triply magic hypernucleus 208+Apb A=2, 8, 

20, 40, and 70. The proton and neutron densities are al¬ 
most not impacted by the hyperons addition, showing a rel¬ 
ative independence of the hyperons with respect to the nu¬ 
cleonic core. The results are in agreement with Ref. 
where the Skyrme parameterisation Sill was used and the AA 
functional neglected. Qualitative agreement, that is stability 
of radial nucleonic mass distribution and non-monotonic evo¬ 
lution of strangeness radial distributions upon increasing the 
number of As, is also obtained with the pioneering work of 
Ref Oil where double and triple magic O and Ca hypemu¬ 
clei have been addressed within a RMF model. This validates 
the present approach and more globally the microscopic pre¬ 
diction of hypemuclei properties. 


C. Bubbles and haloes 

Bubbles and haloes effects have recently been studied in 
A-hypernuclei ll29l . Hyperons are more diffuse in a nucleus 
than nucleons. This could be due to the weaker AA attraction 
compared to the NN one, generating an hyperon saturation 
density about 1/3 smaller than the nucleonic one. This may 
also be due to the AA functional which is much more intense 
in Ref. Il29ll compared to the one used in the present work. 
It emphasizes the importance of taking into account the bond 
energy in order to constrain the AA functional, as depicted in 
section II and III. 

It is well known that the addition of a single-A hyperon 
shrinks the nuclear core ll58l . both from predictions and from 
measurements. It is therefore relevant to study the effect on 
the neutron and proton densities of a large number of hyper¬ 
ons. We find no large effect of the increase of the A number 
on the proton nor on the neutron density in In the case 
of 104+40 A(^j.^ jjjg lambdas act like a glue between protons and 
neutrons and even drive the proton to larger radii, as shown 
by Fig. A similar effect is observed on the neutron den¬ 
sity. This is due to the fact that hyperons attract nucleons at 
larger radii. It should also be noted that the halo effect in 
104+40 A(^j- jg rather small compared to predictions using RMF 
approaches ll29l . 

In the case of bubbles Il2^l59ll^ . there is no strong impact 
of the increase of number of hyperons on the depletion. Fig- 
ure[^displays the proton, neutrons and hyperons densities in 
^'^Si with no hyperons, and the addition of 2 and 8 hyperons. 
As described above, the proton and neutron densities are al¬ 
most not impacted by the addition of hyperons and therefore 
the predicted proton depletion in ^^Si remains. This small in¬ 
terdependence is at variance with relativistic calculations ob¬ 
tained with the RMF approach 1^ . 
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VI. CONCLUSIONS 



r (fm) 

FIG. 15. (Color online) Baryons densities for the magic *‘*"*Cr 40 A 
hypemucleus (solid lines). The neutron and proton densities for the 
**^^Cr nucleus are in dashed lines. 



FIG. 16. (Color online) Baryons densities for hypernucleus 

with A=2 (dot-dashed lines for neutron and proton densities) and 
A=8 (solid lines). The neutron and proton densities for the ^'^Si nu¬ 
cleus are in dashed lines. 


In summary the present results show that the effect of A- 
hyperons on the nucleon core is much weaker than in the RMF 
case. This may be due to the different (NAh-AA) functionals 
that are used in the two approaches as well as to the Fock term 
which is not considered in the RMF case ITTI . It should be 
noted that the present (NA-t AA) functionals are constrained 
with bond energy requirements. 


In this work, we have investigated the relations between A- 
hyperonic data (at low A-density), and high density properties. 
For that we have proceeded in two steps: first we have dis¬ 
cussed the properties of the functionals, based on BHF data, 
and proposed an empirical prescription for the AA term. Then, 
we have explored the hypernuclear chart allowing the func¬ 
tional to vary inside our domain of uncertainty. 

In the first step, we found that the low density part of the 
functional is well determined by the value of the bond en¬ 
ergy, while the high density part is determined by the unique 
free parameter of the model, which controls the density at 
which the AA term changes its sign. As such we show that the 
bond energy controls the low density part of the EoS solely. 
This means that hypernuclear data cannot be used to deter¬ 
mine the high density behavior of hyperonic matter, and other 
constraints from neutron star physics or heavy ion data are 
needed in that region. This general conclusion is certainly 
independent of our model and shall be found also in other ap¬ 
proaches. 

In the second step, the A-hypernuclear chart for even-even- 
even hypernuclei with Z < 120 and A=2, 8, 20, 40, and 70 has 
been calculated using the Hartree-Fock method with Skyrme 
NN functional and designed NA and AA functionals. Six 
(NA + AA) functionals were used, optimized such as to repro¬ 
duce the present experimental constraint of a bond energy of 1 
MeV. The position of the A-hyperdriplines is determined with 
a 7% accuracy, with comparable contribution from the uncer¬ 
tainty on the NN functionals, and on the A related functionals. 
The number of such bound even-even-even A-hypernuclei is 
estimated to 9770 ± 429, leading to an estimation of the total 
number of A-hypernuclei with A <70 and Z < 120 of 491680 
± 59000. 

Significant deviations from Iron-Nickel elements can be 
found for A-hypernuclei with the largest binding energy per 
baryon, especially for A > 20. The spin-orbit weakening of 
the neutron magicity close to the neutron dripline is quenched 
the presence of hyperons. The nucleonic core profile is not 
much affected by the presence of hyperons, allowing for the 
persistence of the proton bubble in ^^Si with additional hyper¬ 
ons. 

The present results shall benefit from the more and more ac¬ 
curate design of the A-based functionals. The measurements 
of A and multi-A hypernuclei, A-density profiles, as well as 
A-E and E-E interactions, would greatly help to provide such 
critical information. 


ACKNOWLEDGMENTS 

The authors thank I. Vidana for fruitful discussions. This 
work has been partially funded by the SN2NS project ANR- 
lO-BLAN-0503 and it has been supported by New-Compstar, 
COST Action MP1304 and the Institut Universitaire de 
France. 




15 


Appendix A: Relation between the bond energy and the 
functional 

For a nucleus hereafter called (A — 1), including one A and 
A — 2 nucleons, its total energy can be approximatively given 
by: 


£(1-1z)«£(^-2z)+ca(A-1), (A1) 

where ca (A — 1) is the single particle energy of the A state 
in the a hypernucleus, taking into account only the NA 
interaction. For the double-A ^Z hypemucleus, we have also 

£(1aZ) « E{^-^Z)+2eA{A) + u[^\pA{A)), (A2) 


where t/^^^(pA(A)) is the A potential term induced by the 
presence of A particles. In the present functional approach, 
we remind that the momentum dependence of the AA inter¬ 
action is neglected, and that is momentum independent. 
Using the local density approximation, it is a function of the 
average A-density in the ^Z hypernucleus, Pa{A)- From the 
BHF framework, it can be shown that EtII 


^ ttW 
£aa=:zPaUI 


(A3) 


Injecting Eqs. (All and (A2i into ( [T6| ), the bond energy 
reads. 


ABaa(A) « 2eA{A - 1) - 2ca(A) - t/f > (pa(A)) (A4) 

It should be noted that, in order to obtain Eq. ( |A4| l, we have 
assumed the independent particle approximation. Since the 
momentum dependence of is neglected, the A-effective 
mass is the same in double-A and single-A hypernuclei, and 
depend only on p^ as in Eq. (34 1 . The difference eA(A — 1) — 


eA{A), which is mostly induced by the rearrangement term in 
the mean field, is calculated to be small. We therefore approx¬ 
imate eA(A — 1) « eA(A) to obtain 


ABAA(A)«-t/f^(pA(A)) = -2^^^^, (A5) 


where Eq. (A31 has been used. 


Kerman and Lipkin assumed that the energy difference (the 
degeneracy raising) between a nucleus and the corresponding 
hypernucleus state where a neutron is replaced by an hyperon, 
is due to i) the mass difference between the A and the neutron 
and ii) the difference between the nucleonic and the hyperonic 
potentials lIMll : 


AE{SAR) = - m{^X) = (Eqa - Vo^) + ('«A - m„) 

(Bl) 

where Vqa and Vo„ are the depth of the hyperon and nucle¬ 
ons potentials, respectively, and A=1 here. 

This corresponds to an excitation energy in the hypernu¬ 
cleus of 


E* = m{^Xl) - m{^XA) = Sn i^X) + {Vqa - Vo«) +Ba (B2) 

where Ba is the binding energy of the hyperon in the single 
A hypernucleus and S„{'^X) the one neutron separation energy 
of the nucleus. 

The questions arises whether the above relation can be gen¬ 
eralised to multi-A nuclei. A straightforward derivation al¬ 
lows to derive the mass (namely the excited state) of the hy¬ 
pemucleus having A hyperons from the mass of the initial nu¬ 
cleus: 


AE{SAR) = A [(VoA - Von) + {ma - m„)] (B3) 

where A > 1. It should be noted that the AA interaction 
is neglected here, as well as rearrangement terms. This could 
result in a significant different form of Eq. ( |B3| l. However, as 
shown below, the energy position of SAR states will mainly 
be impacted by the nucleon vs. A mass difference. 

Eq. ( |B3| l shows that the SAR states in multihyperons are 
expected to display a rather harmonic spectrum, since the last 
term in brackets of the right-hand-side remains rather con¬ 
stant: there is about 30 MeV difference between the neutron 
and the hyperon mean potentials ll26l . and 170 MeV differ¬ 
ence between the neutron and the hyperon masses: 


m{'^Xl) ~ m('^X) -f A.200MeV (B4) 


Appendix B: Strangeness Analog Resonances 

Strangeness analog resonances (SAR) are similar to iso- 
baric analog states, with the transformation of a nucleon into 
an hyperon instead of a transformation of e.g. a neutron into 
a proton. In Ref ||6T1, Kerman and Lipkin studied the SAR 
states between nuclei and excited states of single hyperon hy¬ 
pernuclei. It is interesting to generalise this approach to multi- 
A systems. 


It should be noted that this 200 MeV constant value origi¬ 
nates from the saturation properties in hypernuclei, which in 
turn leads to a constant difference between the nucleon and the 
hyperon potentials. Moreover, in the corresponding hypernu¬ 
clei, these states would correspond to excited states located 
around A. 30 MeV, which becomes unbound already for small 
value of A. In summary, the SAR states have been generalised 
to multi-hyperons hypernuclei. They shall not correspond to 
bound states for most of them, but to resonances embedded in 
the continuum. 
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